Theory of "ferrisuperconductivity" in Ui-^Th^Be 
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We construct a two component Ginzburg-Landau theory with coherent pair motion and incoherent 
quasiparticles for the phase diagram of Ui-^Th^Beia. The two staggered superconducting states 
live at the Brillouin zone center and the zone boundary, and coexist for temperatures T < T C 2 at 
concentrations x c i ~ 0.02 < x < x C 2 ~ 0.04. We predict below T C 2 appearance of a charge density 
wave (CDW) and Be-sublattice distortion. The distortion explains the fiSH relaxation anomaly, and 
Th-impurity mediated scattering of ultrasound to CDW fluctuations explains the attenuation peak. 



PACS numbers: 75.20.Hr, 71.10.Hf, 71.27. +a 

Introduction: Heavy fermion materials continue to re- 
ceive intense experimental and theoretical interest ]jj . In 
particular, the rather spectacular properties of the heavy 
fermion superconductor alloy Ui-^Th^Bcia have defied 
a comprehensive theoretical understanding for over 15 
years. In this system, the temperature (T) vs. concen- 
tration (x) plane displays four distinct ordered phases, 
as shown schematically in Fig. 1 Pressure studies 

, lower critical field data , and the size of the spe- 
cific heat anomalies ||] strongly suggest that all of these 
phases are associated with superconductivity. However, 
below T C 2, local magnetic probe experiments Q reveal an 
anomaly presumed to arise from weak magnetism with 
ordered moment /i < 0.01/is. Longitudinal ultrasound 
attenuation (a s ) data display a peak, two orders of mag- 
nitude above the normal state background, which con- 
trasts with an a s drop in standard superconductors H. 
The normal state of this system is also anomalous. The 
resistivity of UBei3 is reproducibly large (order 100 
cm) at T c and non-Fermi liquid (NFL) like in its temper- 
ature dependence |IJ (a) , j?]] . Fermi liquid behavior and a 
small resistivity are restored by application of magnetic 
field and pressure, suggesting that this residual resistiv- 
ity is intrinsic 0. These results legislate against an in- 
terpretation of the superconductivity of UBei3 in terms 
of pairing within a traditional Fermi liquid quasiparti- 
cle background. Rather, there is electron pair coherence, 
but electron quasiparticle incoherence, similar to the in- 
terlayer tunneling and "stripe" S theories of cuprate su- 
perconductivity. 

In this paper, we propose a new phenomenological 
theory for superconductivity in Ui-xTh^Beia which ex- 
plains many of the observed data and also is consistent 
with the NFL normal state. Motivated by results for 
the two-channel Kondo impurity El (a), and lattice 
pfj| , pl| models, we develop a phenomenology in which lo- 
calized and presumed odd-frequency pairing excitations 
hop coherently on the lattice, but with negative pair hop- 
ping/Josephson coupling from site to site. Our main re- 
sult is that the low temperature phase below T C 2 is a 
ferrisuperconductor, with spatially modulated pair am- 



plitude and coexisting superconducting order parameters 
with zone center and zone boundary center of mass mo- 
mentum. The name derives from obvious analogies to 
ferrimagnetism. Moreover, the ferrisuperconductivity in- 
duces concomitant charge density wave (CDW) and vol- 
ume conserving antiferrodistortive (AFD) order. The 
combined ordering allows us to quantitatively explain the 
anomalies observed in /iSR relaxation 0] and a s [[| ex- 
periments while invoking no magnetic moments in the 
low temperature phase. For a s , our phenomenology is 
constrained by experiment and microscopic theory, leav- 
ing one adjustable parameter to provide an excellent fit 
to both the temperature and frequency dependence. Our 
theory has several testable predictions. 

Microscopic motivations: The two-channel Kondo 
(Anderson) lattice model, which has been rigorously 
solved using dynamical mean field theory in infinite spa- 
tial dimensions Jj], [io|Ji"l| ] , has been shown to display many 
properties consistent with UBei3. This model consists of 
local spin 1/2 moments coupled antiferromagnetically to 
two species of conduction electrons, derived here from the 
quadrupolar Kondo effect JO (a). This model naturally 
explains the NFL thermodynamics and transport prop- 
erties of the normal phase. The calculations also confirm 
the possibility of an odd- frequency superconducting state 
p| , as anticipated from the strong local pairing fluctu- 
ations in the impurity model [0(a), ||. The relevant 
excitations of the normal state are not conventional Lan- 
dau quasiparticles, so coherent single particle transport 
is irrelevant. 

Ginzburg-Landau Theory: To characterize the various 
superconducting states we construct a Ginzburg-Landau 
functional. Assuming the superconductivity is induced 
by two-channel Kondo effect induced pairing interactions 
at the U sites 0, one can argue that the resulting or- 
der parameters should be odd in frequency, and stag- 
gered in real space |p|,pl]-p^| . To identify candidate or- 
der parameters we investigated a negative pair hopping 
model with third nearest neighbor overlaps, and assumed 
a logarithmically divergent local pairing susceptibility at 
the U sites M. The local symmetry of these pairing or- 
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FIG. 1. Sketch of the different superconducting re- 
gions and the corresponding theoretical order parameters for 
Ui-xTh^Beia. Here, x c l ~ 0.02 and x C 2 ~ 0.04. The maxi- 
mum of Tci is T™ ax « 0.6 K. a a peaks at T C 2, and the muon 
relaxation rate increases below T C 2 . In our theory, for T < T C 2 
a CDW and a lattice distortion are generated, which account 
for these data. 

der parameters about the U-sites is A2 Jl5|. Given the 
two-sublattice fee structure of UBei3 , two possible order 
parameters emerge Jlij]: rjr which sits at the Brillouin 
zone center, and the three dimensional order parameter 
r/i which lives at the three X points Qi at the center 
of the square faces of the fee Brillouin zone. These or- 
der parameters transform as the odd parity and even 
parity irreps of the 0\(Fm3c) space group of UBei 3 . 
As shown in Fig. 2, t]y changes in sign from site to site 
on the simple cubic U-sublattice, while r\i has lines of 
alternating phase parallel to Qi. 

The resulting free energy functional invariant under 
time reversal, gauge and space group operations reads: 

3 3 
F = a T \ V r\ 2 + a x <T \r,i\ 2 + P r \vr\ 4 + ft £ h\ 4 

3 3 

+ E h\ 2 h\ 2 + E vhf + Ft-x , (i) 

3 3 

Fr-x = € E WlTkl* + CE^f + c-c) ■ (2) 

i=l i=l 

As usual, we assume linear temperature dependence in 
the quadratic coefficients, viz., a x = a\(T — T c x) and 
a T = a 2 (T - T c r). 

Our conjectured phase diagram is shown in Fig. 1. For 
pure UBei3 , analysis of observed anisotropics in the up- 
per critical field data [[16[ lead us to identify the ^-phases 
as likely candidates ]17| . Pressure experiments strongly 
suggest that for x c \ < x < x C 2 the lower temperature 
phase represents a continuation of the low concentration 
line Next, we conjecture that Th doping produces a 
relative stabilization of the rjr phase at intermediate con- 
centrations, so that r\i and rjr coexist for x c \ < x < x C 2 



and T < T C 2- Given X and F point order parameter 
coexistence, we call this state ferrisuperconductivity in 
analogy with ferrimagnetism. 

The parameters, (if , £ and £ in Eqs. ([!]) and (|2|) are 
not known a priori, but they are largely constrained by 
the physical properties of Ui-^Th^Beia. Since there is 
no CDW formation, time reversal symmetry breaking, or 
splitting of the transition under applied stress in pure 
UBei 3 |5), the A X ' S must satisfy > [3 X and (3g < 
~ Pii yielding a phase with a single r)f condensed. 

The nature of the coexistence phase depends on the 
sign of £■ For £ > we find a time reversal broken co- 
existence phase, where microscopic current loops induce 
interstitial antiferromagnetic order, a possibility noted 
previously for UPt3 Jl8) . This could explain the observed 
/iSR anomaly. However, it is more difficult to explain the 
observed large peak in a s at T c2 || . For ( < the mixed 
superconducting state is completely different: in this case 
a CDW and a staggered lattice distortion appear below 
T C 2 and the cubic symmetry of the lattice is lowered to 
tetragonal (see Fig. 1). The CDW arises because the 
phasing of the rjr , r\i order parameters corrugate the pair 
amplitude locally, producing a ±|r7r||r7i| alternation from 
plane to plane perpendicular to Qi (see Fig. 2). 

The anomalies of the second transition can be ex- 
plained consistently and easily for Q < 0: the lattice 
distortion can produce the observed increase in the /zSR 
relaxation, while the peak in a s results from coupling to 
large induced CDW fluctuations. 

Henceforth, we concentrate on the coexistence phase of 
Ui-^Th^Beis. To estimate the amplitude of the CDW 
(gi) and the displacement of the face Be atoms (ui) we 
first construct the corresponding symmetry allowed cou- 
pling term in the Ginzburg-Landau functional: 

3 2 3 2 

^g-rj-u—Xg 2_^i 2 ' J 2 

i=l i=l 
3 3 

+ QiUi + Qi(rirri* + c.c.) , (3) 

i=l i=l 

where the Uj's and the Qi's form three dimensional X^ = 
T2 ® X^ X point reps. The coupling 7 w 0.4 eV/A and 
the spring constant K u ~ 2 10 4 K/A 2 can be obtained 
from microscopic Thomas-Fermi theory calculations and 
from the bulk modulus, respectively, while A sa 10 K is 
estimated from the relative shift of T c2 and x" 1 - 900 K 
is extracted from fits to the experimental data. 

As we discuss below, from the ^(SR experiments we 
estimate a T — OK lattice displacement amplitude uq — 
0.4 A for Be on planes perpendicular to Qi. Using Eq. (|^) 
a simple calculation yields an associated CDW amplitude 
of Qi < 1.5, where we simplistically use uq in the spring 
energy to produce an upper bound. This CDW ampli- 
tude is quite large at first sight. However, the charge is 
distributed among the 24 face Be sites surrounding the 
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FIG. 2. Schematic depiction of proposed ordering in 
Ui-xTh^Beia below T C 2- Large circles are U atoms, small 
circles represent face Be atoms, viewed on edge. Q\ repre- 
sents the propagation direction for the assumed X point pair, 
charge, and displacement order. The relative sign of the CDW 
order is shown by ± notation on the U atoms. Arrows repre- 
sent associated displacements of the Be atoms. 

U ion, resulting in an excess charge of ~ 0.1 e/Be. uq 
is somewhat large, but in view of the large lattice con- 
stant a » 10 A, is in reasonable agreement with observed 
displacements in CDW systems Q. Moreover, it is stag- 
gered and volume conserving (to first order). 

The displacement does result in an increased fiSK de- 
phasing rate below T C 2, with reasoning as follows: The 
muon sits in the center of the Be(II) plaquettes (20). Be 
plaquettes on faces parallel to Qi are alternately com- 
pressed and expanded. The dominant spin-dephasing 
mechanism of the muons is due to the dipolar coupling 
to Be nuclei, and varies in quadrature as R~ 6 , where 
R is the muon-nearest-Be separation. Thus, muons 
in expanded (compressed) plaquettes will experience re- 
duced(enhanced) relaxation. Averaging over inequiva- 
lent stopping sites, leaves only even powers of itj, and 
coupling to compressed plaquettes dominates leading to 
an increase in dephasing. The resulting change in the 
muon relaxation rate out is given, to leading order in 
uo, by 

*Z™ « B ul = B - \vrW ■ (4) 

with B « 0.4A 2 . Taking uq = 0.5 A yields the estimated 
zero temperature relaxation rate increase of ~ 10%; ex- 
tending Eq. (4) to all orders in uq gives uo ~ 0.4A 
= 0.04a, in rough agreement with NMR imposed con- 
straints plj . The leading order scaling with 77? is in ex- 
cellent agreement with the experimental observations . 

Now we analyze the anomalous peak in a s . As shown 
in Ref. [p|, a s is practically the same for (100) and (111) 
propagation directions, ruling out the domain wall dis- 
sipation mechanism of Refs. p^p3| . The original exper- 
imental reference II (and recent interpretation of ther- 
mal expansion data || ) suggest that the transition below 
T C 2 is purely magnetic, citing large moment Gd and Mn 
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FIG. 3. Fig. a: Fit to the experimentally observed ultra- 
sound attenuation peak. To evaluate the integral in Eq. ([?|) 
we used a cutoff E c ~ 0.3 K. The transition temperatures are 
Tci = 0.47 K and T c2 = 0.37 K. Fig. b: Leading correction 
to the ultrasound attenuation. The wavy lines indicate the 
acoustic phonon propagators, while the dashed line with the 
cross stands for impurity scattering at the Th sites. Bubbles 
represent CDW and pairing susceptibilities. 

based metallic and insulating magnets with a s ^ mag wlO 
dB//isec at these frequencies. However, for magnetic 
coupling, the two primary contributions to attenuation 
are from energy relaxation (which varies like /1 2 ), and 
order parameter fluctuations (which varies as /1 3 ) p4[ . 
These mechanisms are incorrect in temperature depen- 
dence and magnitude: given fx <0.01/x_b from /iSR [||, 
we find a s < (.01^s/.5^s) 2 a S!mas = 4x 10~ 5 dB//xsec, 
four orders of magnitude too small. 

In view of these results, and the above interpretation 
of the /zSR anomaly, we propose that the a s peak is due 
to the strong interaction of the impressed phonons with 
CDW fluctuations at the Th sites. Sound waves propa- 
gate with a momentum q sa 0, and therefore can linearly 
couple to the CDW fluctuations at the X points of the 
Brillouin zone only via the Th impurities at sites {Rimp}- 
The strongest coupling mechanism we find results from 
the destruction of the superconductivity and thus the 
CDW at the non-magnetic Th sites, with coupling 
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(5) 



where u denotes the acoustic displacement field and g ~ 
0.8 eV is estimated from microscopic model calculations. 
a s is obtained from the imaginary part of the acoustic 
phonon self energy, shown to leading orxder in Fig. 3. 
Replacing the CDW bubbles by the renormalized static 
susceptibility, \g = Xe/0- ~ l 2 Xg/K)i an d assuming a 
quadratic dispersion of the form w(k) = C\T — T c ?\ + 
D\k\ + D\\k 2 with a cutoff energy E c ~ IK we can 
integrate out the internal momenta and we arrive at the 
following expression: 



a s w Im{n(o;)} = Cx g 



2 X 2 \vr\ 2 ^ 



I = - 



dx 



Mc 2 r 
4T 2 lu 2 



(6) 



(7) 



/0 \\a x \ + E c x) 2 +AT 2 u 2 1 

with C a number of the order of unity. To derive Eq. (]?]) 
we assumed a hydrodynamic damping 1/r « ujT for 
the superconducting fluctuations. T is the only unde- 
termined parameter in Eq. (0); all others are fixed by 
microscopic model estimates or experimental constraints 
p5[ . Substituting the different couplings into Eq. (||) we 
estimate the prefactor to be 20 r _1 dB//is, which is in 
very good agreement with the experiments (away from 
criticality, the overall magnitude is considerably reduced 
by the integral in Eq. (6)). Note that: (i) a s is asym- 
metric about the transition due to the ij 2 factor; (ii) for 
lu -> 0, a s ~ |T— T c2 |~ 1/2 away from T c2 ; (iii) for T = T c2 
and oj — > 0, a s ~ w 3 / 2 . We obtain a fit of remarkable 
quality to the a s data of Ref. f|, as shown in Fig. 3. 

Predictions and further considerations: Our theory ad- 
mits a number of falsifiable predictions below T c %: (1) 
The Be displacements with magnitude < 0.4A should be 
easily observable in neutron scattering experiments. (2) 
The impurity scattering mechanism produces an isotropic 
angular dependence for a s above T c2 and zero attenua- 
tion for transverse ultrasound, consistent with the rough 
agreement between longitudinal attenuation coefficients 
along the (100) and (111) directions measured for differ- 
ent samples of comparable nominal composition Q . We 
anticipate an order of magnitude difference between lon- 
gitudinal and transverse attenuation. (3) The a s peak 
will be maximized at the concentration x f=a 0.03 and 
drop significantly in magnitude near x c i,x C 2- Assum- 
ing insignificant variation of electronic parameters as a 
function of x, then a s (a;)/(T c i(a;) — T c2 (x)) should vary 
as x according to Eq. (6). This may provide a test of 
the impurity scattering mechanism. (4). The observed 
NMR line broadening below T c2 of Ref. |2l| was origi- 
nally attributed to the flux lattice, but this is ruled out 

The broad- 



by subsequent penetration depth data |26|. 
ening is roughly consistent with field gradient shifts due 
to our proposed Be sublattice distortion. 
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